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Abstract 



Production of lepton pairs in the 7p e~ e~^p reaction is studied at photon 
energies up to 1 GeV. We show that even if lepton charge is not measured there 
are extensive kinematical regimes where the nuclear process dominates over 
the Bethe-Heitler contribution. The decomposition is performed of the unpo- 
larized cross section for the virtual Compton scattering in terms of response 
functions. These are expressed in terms of a polarization density matrix of 
the virtual photon and are shown to be sensitive to the properties of baryon 
resonances in the first and the second resonance regions. In the analysis a 
unitary K-matrix model is used, based on an effective Lagrangian including 
nucleon, mesons, and baryon resonances with masses up to 1.7 GeV. Results 
of the model are compared with data for real-photon Compton scattering and 
predictions are given for observables in photoproduction of lepton pairs. 

1996 PACS numbers: 13.30.-a, 13.40.-f, 13.40. Gp, 13.40.Hq, 13.60.Fz 

Key Words Virtual Compton scattering; Real Compton scattering; Baryon resonances; 

K-matrix approach 



I. INTRODUCTION 



Compton scattering (CS) on the proton with real and virtual photons is a fundamental 
process which provides information on the internal structure of the proton and its excited 
states. Virtual CS where the initial photon is space-like, i.e. the reaction e~p e^jrf , has 
recently attracted considerable attention. Below the pion production threshold it allows, via 
an interference of the Bethe-Heitler (BH) process with the proton virtual CS, to measure 
the generalized polarizabilities of the proton. For a review we refer to . 

Virtual CS in the time-like region, the reaction e~e'^p, is of interest for com- 

plementarity reasons. The process is sensitive to the nucleon electromagnetic (e.m.) form 
factors in the region < < 4m^ (where m is the nucleon mass, k is the photon momen- 
tum), which are not accessible in elastic electron-proton scattering or e~^e~ annihilation to 
pp. This possibility was explored in Ref. 0. In the present paper the main interest is to 
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investigate dependence of the polarization density matrix of the virtual photon on various 
contributions to the nuclear matrix element. 

In certain parts of phase space the cross section is dominated by the BH process. Only 
when lepton charges are measured this can be turned to an advantage by measuring an 
e^e~ asymmetry which is directly proportional to the virtual CS - BH interference. This 
was suggested long ago @] and more recently elaborated in 

In this paper we study the situation where the electron and positron in the 7p — >■ e~e^p 
reaction are not distinguished which is the case for experiments which are being done at 
MAMI (Mainz) [Q. Under such conditions, because of the different charge-conjugation 
parity of virtual CS and BH amplitudes ||^, the above interference vanishes and one is 
left with the incoherent sum of virtual CS and BH contributions. We show that there are 
large kinematical regimes where the more interesting virtual CS on the proton is dominant, 
mainly at backward angles for the virtual photon. In these conditions the (unpolarized) 
7iV — i> e~e~^N cross section can be decomposed in terms of response functions (RF)s, like 
is usually done for exclusive electron scattering on nuclei (see, e.g. 0]). In addition to 



the transverse RF, which can also be measured in unpolarized CS with real photons, there 
are three more RFs for dilepton production. These are directly related to polarization 
properties of the virtual photon. To make this link more transparent we apply to the 
jN — »• 7*iV reaction the density matrix formalism developed |11| for the photoproduction 
of vector mesons on the proton. 

Differential cross sections and RFs for the 7J9 e~e^p reaction are calculated in the 
unitary K-matrix model previously developed ||12|J13| for pion-nucleon scattering, pion pho- 
toproduction and real-photon CS. Before applying the model to virtual CS we compare its 
predictions to data for real CS in the energy region up to 1 GeV. 

In the calculation of virtual CS we concentrate on the 1st and the 2nd resonance regions 
which are clearly seen in the cross section for real CS. The cross section and elements 
of the polarization density matrix pw are calculated for two energies corresponding to 
excitation of the P33(1232) and the Di3(1520) resonances. We study effects of different 
reaction mechanisms on these observables. In particular, for the A region it is of interest to 
investigate the effect of t-channel exchange of the a meson. 

At photon energies of about 700 MeV photon invariant masses up to ~ 540 MeV become 
accessible. In this energy regime several contributions are important, in particular, the 
Di3(1520) and 5*11(1535) resonances, and a exchange. Due to their different influence on 
the density matrix elements pxy these contributions can be distinguished in an experiment 
provided RFs are extracted. Cross section and other observables are also sensitive to the 
Gs{k'^) coupling in the 7A^iI'i3(1520) vertex which is specific for virtual photon and does not 
contribute to real CS. 

The paper is organized as follows. In sect. || we specify variables and give the expression 
for the exclusive cross section. The decomposition of the virtual CS cross section in terms of 
RFs is presented. Relations between RFs and polarization density matrix elements are es- 
tablished and the limit of "almost" real photons is considered. We briefly discuss ingredients 
of the dynamical model used and compare calculated cross sections and polarization observ- 
ables with data for real CS. In sect. |T| kinematical regimes are investigated where the BH 
contribution can be neglected compared to virtual CS. Results for response Ws = Wt + Wl 
and elements of the photon density matrix are presented. Conclusions are given in sect. 
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In Appendix A expressions for RFs for virtual CS and polarization observables for real CS 
in terms of helicity amplitudes are collected. 



II. FORMALISM 
A. Cross section for e~e^p reaction 

— * 

If g = {qo= \<i\-,<i) is the momentum of the incoming real photon and k = {ko, k) = k_+k^ 
the momentum of the outgoing virtual photon (see Fig. (|l|)), the fully exclusive cross section 
in the cm. frame can be written as Q 

da{e-e+) _ ^ , . . 



K^e-e+\^BH + ^VCsl , (1) 



dM^dn^dni 

where My = y/W is the photon invariant mass and /Ce-e+ is the kinematical factor 

^ m^ml \k\ 4\lf 

2(27r)5s \q\ M^(3^ko ' ^ ^ 

Here m{me) is the proton (electron) mass, (3 = {1 - Aml/M^f'^ is the e~ (e+) velocity 
in the virtual-photon rest frame, and s is the invariant energy squared. The average over 
initial and sum over final polarizations of all particles is implied in Eq. (|l]) as well as in the 
following. 

To describe the e~e^ pair the relative 4-momentum / = \{k- — A;+) has been intro- 
duced [|15|. The orientation of / is determined by the polar angle di and the azimuthal (or 
out-of-plane) angle (f)i which are defined in the frame with OZ axis along k and OY axis 
along k X q. The magnitude of I is given by 

with ko = [s — m'^ + M^)/2yG and k^ = kl — M^. The sohd angle differential dfli stands 

for d cos 6id(j)i, and dQ^ = dcos6^d(f)^, where 6^ is the angle between vectors k and q. 
The azimuthal angle 0^ is superfluous and can be chosen zero. For more details about 



kinematics we refer to previous papers on e~e production in capture reactions [0,0] and 
NN bremsstrahlung ||17|| . 

The Bethe-Heitler (Abh) and the virtual CS [Ayes) amplitudes are functions of five 
variables: s, M^, 6^, 6i and (pi. Before proceeding further one can exploit the different charge- 
conjugation properties of these amplitudes |^. Interchanging variables describing e~ and 
e"*", which means I — > — / and interchange e~ and e"*" helicities, leads to the relations 



^We follow notation of Ref. 
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\ABHiOi,(f)i)\^ = \Abh{ti-0i,7i+''^'^ 



D 

\AvcsiOi,(l)i)\^ = \AvcsiT^-ei,n + (f)i)\^ 
A*j^HiGi, <Pi)Avcs{Oi, (pi) = ~A*sHiT^ -OuT^ + (t)i)Avcs{T^ ~Oi,n + 0;) , (4) 

where 6i varies from 0° to 180° and only the dilepton angles are indicated. Therefore, if 
the electron and positron cannot be distinguished in an experiment, the cross section is the 
incoherent sum da{e~e~^) = da{e~e~^) + da{e~^e~) , which results in 

^^^"""^^ -2IC^-^.{\AnH\' + \Aycsn, (5) 



dM^dVt^dVti 



where the interference term drops out. 

Another interesting observable, the e~e^ asymmetry, which is proportional to the virtual 
CS - BH interference (see @-0)7 will not be addressed here. Measuring this would require 
detection of the charges of the leptons. 

The expression for Abh can be read off the diagrams g,h in Fig. (|l|). After some algebra 
it takes the form 

Abu = Q-Ak-)M^^ -ir-q^ + 2r-q " 2^77^ ^^^'-^^ 

xuip',A')[F,iQ')Y + ^—^F,iQ')]uip,A)eM , (6) 

2m 

where u{p, A) {u{p' , A')) is the spinor of the initial (final) nucleon with momentum p {p') and 
helicity A (A'). In Eq. (P) e is the proton charge, u{kJ) and f (fc+) are the lepton spinors 
(omitting spin indices), and e(A) is the photon polarization vector for helicity A. The nucleon 
e.m. form factors i^i,2(Q^) depend on the momentum Q = p' — p. 
The virtual CS matrix element is 

Aycs = ^u{k_)^AK)J\ (7) 

where the e.m. current is expressed through the CS tensor M^^ 

J'^ = eu{p\ A')M^''u{p, A)e^(A) . (8) 

Some of the ingredients in the calculation of M^'^ are discussed in sect. P C| . 

Note that the real CS matrix element {Abcs) can be written in this notation as 

Arcs = e\';{X')u{p',A')M'^''u{p,A)e,{X) (9) 

with e'(A') being the transverse polarization vector of the final photon which is real in this 
case. 



B. Response functions and polarization density matrix 

In kinematics where virtual CS is the dominant process it is of interest to decompose 
the cross section in terms of response functions (RF)s, similar to what is common practice 
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in the space-hke region. We will not go in much detail because the derivation is similar to 
e~e"'" production in other reactions ||T5|-17 . 

After summing over lepton polarizations and making use of gauge invariance k ■ J = 
one obtains in the cm. frame 



do-(e e"*") 
dM^dVt^dVti 



2S 



S 



2TT^sM^f3^ko \q\ 
WtXt + WlXl + WttXtt cos 



+ WltXlt C0S( 



(10) 
(11) 



where a is the fine-structure constant, Wi = Wi{s,9.y, M^) are RFs for i = T,L,TT and 



LT , and the factor 2 in front of S comes from the definition Eq. 
polar angle 6i is contained in the factors Xj 



). Dependence on the 



XT 



XtT 



sin^ Oi 



XL 



P 



sin^ Oi 



xlt 



M^ko 



sin 201 



(12) 



where |Zj, which also depends on 9i, is defined in Eq. (^) 



RFs are defined in terms of the space components J* of the e.m. current as 
1 



T 



Wr 



polar. 



1 k'^ 

___ I ? |2 

Akl V 

u polar. 



W' 



TT 



E \Jy 



J T. 



LT 



polar. 



1 

polar. 



(13) 



where the sum runs over the proton and initial photon polarizations. 

It is interesting to note that certain combinations of the e.m. current vanish due to 
reflection symmetry {OY —OY) with respect to the scattering plane OXZ: 



E '^{JiJy)= E ^{JiJy)= E '^{J:Jy)= E ^{J:Jy) = ^- 

polar. polar. polar. polar. 



(14) 



These constraints are specific for a two-body final state such as iV -|- 7* and may be useful 
for checking consistency of a model. 

Experimental separation of RFs can be performed using the different dependencies of 
the kinematical factors in Eq. (0) on the dilepton angles. In Refs. |18] an example of such 
a separation for the case of A^A^ virtual bremsstrahlung is discussed. 

To obtain some further insight we relate the above RFs to elements of the polarization 
density matrix p of the virtual photon. The formalism of the polarization density matrix 
has been used before in production of vector mesons p.uj (see Refs. |lT9| , p!Tlj20| , |2l|| ). One 
introduces the 3x3 density matrix for a spin-1 particle 



/ Pll 

\pi-i 



Pio Pl-l 
Poo — Pio 
-Pio Pll 



(15) 
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in the helicity basis A, A' = 0, ±1. It satisfies the conditions of hermiticity paa' = Pa'Aj parity 
conservation pxy = {—l)^~^p_x_\', and has a unit trace (pn = |(1 — poo))- Poo and pi_i 
are thus real quantities while the element pio is, in general, complex. 

The form Eq. (|15D is valid only for a reaction with a two-body final state (like + 7*), 
in which the initial hadron and photon are not polarized and the polarization of the final 
hadron is not measured. If the initial photon in — 7*A^ is (fully or partially) polarized 
the general density matrix is more complicated |TI|. 

From the density matrix and the 76" e+ vertex, which is well-known from QED, one 
can obtain the angular distribution F{6i,(f)i) of 7* decay into e~e^. In the virtual-photon 
rest frame it takes the form (see Ref. |2^, ch. 4-13 [| ) 

F{ei, 0,) = ^ { 1 - [ ^(3poo - 1)(3 cos^ 01 - 1) 

— 3pi_i sin^ 61 cos 2(f)i — 3\/23?pio sin 26i cos (pi]} , (16) 

where P is defined after Eq. (|) and {9i, (pi) denote the dilepton angles. 

Expression Eq. (pUD in terms of RFs, boosted to the 7* rest frame, can be cast in the 



form (details on the angle transformation can be found in Refs. ||T^,[T7 



S = Wt{1 - sin^ Oi) + Wl{1 - cos' 9i) + Wtt\p^ sin' Oi cos 

+ W^LT^/?^ sin 2^; COS 0; 

= (l--)iy,{l-^[ (3 COS ^^-l)-^sm ^,cos20, 

where Ws = Wt + Wl . It is seen that the angle dependent part of S has the same structure 
as F{9i, (pi) in Eq. (|16|). Comparing terms with the same angle factors one obtains 



Wl Wtt ^ Wlt 

P''=w-s^ P'-' = 2Ws^ ^P'' = mrs- ^'^^ 

Note that the imaginary part of pio does not enter in the 7* —>■ e~e~^ angular distribution 
(cf. [|ll| for 7A^ — > pN —>■ TCTcN). For completeness we give an expression for Qpio through 
an additional (independent) function Wlt- 

^pio = ^ , Wlt = E 2v^^ {.rM . (19) 

^ ^ polar. 

Wlt and Wlt can be regarded as the real and the imaginary parts of one complex function. 
Thus only three (out of four independent) components of the photon-polarization density 



^There is a misprint in Eq.(13.14) of this reference: in front of oi there should be an additional 
factor of 1/2 
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matrix in addition to the cross section can be determined from measuring Wt, Wl, Wtt, 
and Wlt in the — e~e~^N reaction with unpolarized leptons. 

Since |Wrr| < Wt < Ws (see definitions in Eqs. (0)), the following inequalities are in 
order 



< Poo < 1 



, 1. X 1 

|pi-i| < -(1 - Poo) < 2 • 



The upper bound for pio was obtained in Ref. [0 using the Schwartz inequality 

1 / — : r 1 / — : r 1 

2pi-i) < -^V^oo(l - Poo) < • 



IPioI < 2 VPoo(l 



Poo 



2V2 



(20) 



(21) 



Additional constraints have recently been derived in . 

Sometimes (see [^,^) the density matrix for the spin-1 particle is expressed in terms 
of tensor and vector polarizations 



V 



1_+ ^t20 
2^11 ~ V 2^21 



2^11 



2^21 



22 



1 - V2t 



20 



-tu + 



22 



fill + ^J^hl - I ^ 



2*- 

1 + ^^20 



21 



(22) 



where the tensor polarizations ^20? ^21 and t22 are real, and tn is purely imaginary and related 
to the vector polarization py via tn = —i^py. Relations between these polarizations and 
elements pw follow from Eq. (|T5|) and Eq. 



Poo 



5^Pio = —-^hi 



(23) 



Using Eqs. ([ISD and Eqs. (123|) one can express the polarizations in terms of RFs 

1 Wt-2Wl V3Wtt 



t 



20 



t 



21 



V2 Ws 
VS Wlt 
2V2 Ws 



22 



t 



11 



2 Ws 
: VS Wlt 
'2V2 Ws 



(24) 



As a last point it is interesting to consider the limiting case of small invariant masses 
My close to the threshold value M*'' = 2me. In this case the longitudinal component of the 
e.m. current becomes negligibly small, and elements poo and pio vanish. As follows from 
Eqs. dH), t2o I/V2, ^21, tu 0. In this limit Eq. (||) reduces to 



1 



Wtt n 
Wt 



^ Wt 




1 



(25) 



From a physical point of view this limit is close to real CS. A real photon is, in general, 
described by the 2x2 density matrix in the helicity basis [O] 



P7 



1 / 1 ± /c -h exp{-2t^) 
2l-/Lexp(2#) iTic 



(26) 
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where Il{Ic) is the degree of hnear (circular) polarization, ip is the angle between the 
direction of linear polarization and the OX axis, and + (— ) stands for the right (left) 
circular polarization. Comparing Eq. (|26| ) with Eq. ( pSj) (dropping there the superfluous 
line and column with elements equal to zero) one finds 

\Wtt\ 

_ f , if sign(Vrrr) = -1 , 
^ ~ I if sign{WTT) = +1 ■ ^ ^ 



2 ' 



Thus the "quasi-real" photon has no circular polarization, and may have a linear polariza- 
tion along the OX axis (in the scattering plane), or along the OY axis (perpendicular to 
the scattering plane) depending on the sign of the interference Wtt- Note that this is a 
particular feature of any reaction with unpolarized particles in coplanar kinematics and can 
be understood based on reflection symmetry with respect to the scattering plane. 

The photon asymmetry ||2^, as can be measured in real CS with a linearly-polarized 



photon beam, is related to and W^j. (the superscript "r" denotes taking the real CS 
limit). If there are no bound states in the energy regime of interest, the RFs for e~e^ 
production behave smoothly as functions of Q and the threshold values Wt, tt at = 
2me can be used to calculate rpj,. Making use of Eq. (^61) and time-reversal invariance 
one obtains for the real CS cross section 

^ = ^ (^^ - ^TTh COS 2^) , (28) 

where the angle ip specifies the direction of the incoming photon polarization and all other 
particles are not polarized. The cross section for unpolarized photons is expressed via Wj 



rr 

and the asymmetry is proportional to the interference W^j.: 



rfcr_L(7) + c?o-||(7) Wr 



T 



Calculations of the cross section and RFs for virtual and real CS are performed in the 
helicity formalism (see Appendix A where details are presented). 

It is of interest to make a link between the representation in Eq. ( p!TD and the longitudinal- 
transverse (L-T) decomposition established for electron scattering on nuclei. At invariant 
masses considerably larger than 2me one can rewrite Eq. ( p!TD as 

S = 2xt[ {^Wt) - e+ {-Wl) - e+ cos 20^ {-]^Wtt) 



+ 77^2£+(l-£+)cos0, (i^Wlt) ] , (30) 



^In principle, kinematical conditions exist when the electron interacts strongly with the proton 
via the Coulomb attraction and may be captured in a hydrogen atom. However, this happens in 
a tiny region of the phase space and is neglected in the present model. Likewise the Coulomb 
interaction between the e~ and the e"^ and possible formation of positronium is not considered 
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where t] = sign (sin 25^), and e+ = (1 — xt)/xt which satisfies the conditions < e+ < 1. 
Expression Eq. ( ^0|) is similar to the L—T decomposition used in electron scattering (see e.g., 
p!0| , |25[] ) To obtain exactly the same form one can make in Eq. ( pO] ) the formal substitutions 
k"^ — >• — |/c^| and — >• — This results in 



S = 2xt[Rt + £- Rl + £- cos 2(1)1 Rtt + V 2^- (1 + £-) cos i?LT ] , (31) 
where RFs for the electron scattering, denoted by are related to Wi via 

Rt = \wt, Rl = -Wl, Rtt = ~Wtt, Rlt = -^Wlt. (32) 

Comparison of Eq. ( |3TD with the L — T decomposition in (Eqs. (22,23)) suggests that e_ 
has a meaning of degree of transverse polarization of the space-like photon. The interference 
Rlt is of course a real function in this case since the factor i merely accounts for y/k^ — > 



iylk"^] in the definition of Eqs. ([TB|). 

C. Description of the model and comparison with data for real Compton scattering 

As mentioned in sect. | we employ the K-matrix model which is unitary in the chan- 
nel space vrA^ © 7A^. The tree-level K-matrix includes the four-star baryon resonances: 
Pii(N1440), Di3(N1520), Sii(N1535), Sii(N1650), P33(A1232), S3i(Al620) and D33(A1700) 
in the s- and u-type diagrams. The t-channel contributions come from the exchange of a 
and p mesons (for 7iN scattering), vr, p, and u mesons (for pion photoproduction), and vr", 
7] and a mesons (for CS). The amplitude for CS is shown in Fig. (|l|), diagrams a - f . For 



more details about the model we refer to Ref. |jT3|. In order to calculate the virtual CS the 
longitudinally polarized photons have also been taken into account. 

In this paper we concentrate on the P33(1232) and the Di3(1520) resonance regions. The 
^NR vertex for the e.m. couplings of the spin-3/2 resonances is chosen as follows [^^ -p9 



Vn.^r'^ = ^ [G,{ky^p{z,)js - '^^e^p{z^)ps - ^^e^p{z,)h]T {k^e' - k't^) , (33) 

and yRQ^jv7 = — 7o^Ar7^Ra7o • In EQ- O) P [k) is the nucleon (photon) momentum and 
e is the photon polarization vector. For isospin-1/2 resonances Gi{k'^) = [gi^\l + t^)/2 + 
(^•"•'(l — T^) /2]FvMD{k'^) 5 while for isospin-3/2 case Gi{k'^) = giT^FvMD{k'^) , where T3 is 
1/2 3/2 isospin-transition operator and FyMDik"^) is the form factor in the extended 
VMD model ||30|. Note that the G^lk"^) term contributes only for virtual photons. Further, 
F = (75,1) for the resonances (3/2+, 3/2^), and the tensor Oapi^Zi) = g^p + ailalp with 
Oj = —(1/2 + Zi) contains the off-shell parameter zi p6| , p7| , where i = 1,2,3. In the 



present work we show calculations for two parameter sets "A" and "B": A) = for all 
spin-3/2 resonances, B) ai = 0.12, 02 = 0.5, as = for the Di3(1520), and = for other 
resonances. The fitted values of the couplings {gi, g2) are: (4.9, 5.27) for the P33(1232), 
(7.25, 7.9) for the ^13(1520) (on the proton), and (1.74, 4.75) for the ^33(1700). 

The cross section for real CS on the proton and polarization observables are presented 
in Fig. (^ for photon energies up to 1 GeV. The cross section in the A-resonance region 
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agrees with the data as well as with other calculations which are based on dispersion re- 



lations |3T|j3^j3J], and an effective Lagrangian For comparison we show results of the 
calculations from Ref. |^T| (dotted lines) as they extend to high energies. In Ref. the 
imaginary part of the Compton amplitude is calculated directly from pion-photoproduction 
data and the real part is calculated using fixed-t unsubtracted dispersion relations. Our 
model and |3ll] give the same results for the photon asymmetry at the A resonance, 
however they predict a different slope. Similar disagreement was noticed in — it may 
point to the importance of analiticity constraints in the calculation of the Compton ampli- 



tude. The data for the photon asymmetry at 6^ = 90° and energy El^ > 500 MeV 



are not reproduced by the calculation with the set "A" . The description improves with the 
set "B" (dashed lines) though both parameter sets give very close results for cross sections 
(top panels). For the photon angle 120° the difference between solid and dashed lines for 
diminishes and the data do not allow to distinguish between the two parameter sets. 
The recoil-proton polarization P^^ (Fig. (0), bottom) vanishes below the pion-production 
threshold since the imaginary part of the amplitude due to vrA^ rescattering is zero Q. The 



behavior of P^^ is in agreement with the present calculation as well as with the model pT 
although the error bars on the data are large. 

Above ~ 900 MeV the present model loses predictive power as tails of the higher reso- 
nances (not included in the model) will start to influence the cross section. 

III. RESULTS OF CALCULATIONS AND DISCUSSION 

A. Cross sections for virtual Compton scattering 

Calculations of the exclusive differential cross sections for 7p — e~e~^p at two photon 
energies are shown in Fig. (^. The energy E^'''' = 320 MeV corresponds to the A resonance, 
and Ei^^ = 700 MeV roughly corresponds to the region of the Di3(1520). At fixed energy 
the cross section depends on four essential variables. In presenting the results we will 
therefore fix some of them. The virtual-photon scattering angle 6^ is chosen in the backward 
hemisphere as our calculations indicate that this provides favorable conditions to suppress 
the BH contribution. 

In Fig. (^) the azimuthal lepton angle was fixed in plane, i.e. (pi = (coplanar kinematics). 
For the case when e~ and e"*" are not distinguished the angle combination {(pi = 0, 6i} is 
equivalent to {(pi = it, it — 6i}. The results are shown as function of the polar angle 6i at 
several values of the photon invariant mass. 

The feature of the BH cross section (dashed lines) in coplanar kinematics is a pronounced 
peak which develops when either electron or positron moves along the direction of the initial 
photon. The propagator of the lepton in the BH amplitude in Fig. (^ (diagrams g, h) 
becomes large in this case. As is seen from Fig. (^, with increasing the BH peak is 
shifted towards larger angles 9i, and for = 500 MeV (Fig. (^, right, bottom) it occurs at 



^Contribution to the imaginary part from the rescattering is of the order and thus negligibly 
small 
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the angle 6bh = 32°. The maximal photon invariant mass at a given energy is determined 
from the relation M^°^ = ^/s — m and corresponds to production of leptons with energies 

e_ = e+ = ^M^"^^ and momenta k- = —k^. At this so-called kinematical limit the position 
of the peak reaches the maximal angle 6^ff =11 — 9^ which, for instance, for E}^^ = 700 
MeV and 9^ = 135° takes the value 45° at M^'^^ = 543 MeV. 

An example of a calculation in non-coplanar kinematics is shown in Fig. (^. It is seen 
that the BH peak disappears as the momentum of e~ or e"*" cannot be coUinear to the photon 
momentum in this case. The (pi dependence of the nuclear CS cross section is governed by 
the interference RFs which stand in front of cos (pi and cos 20; in Eq. (pH]) . The (pi dependence 
of the BH contribution is more complicated. 

At small invariant masses the nuclear CS cross section (Fig. (|^), solid lines) shows a 
distinct peaking at forward and backward angles, while at large invariant masses it tends 
to flatten. This behavior can be understood from Eq. ([T0|) . Assuming that (M^/2me)^ ^ 1 
one can show that the relative e~e~^ momentum obeys the condition 

1 ^ r |M^(sin2^, + M2/fc2)-i/2^ if M^«|fc|, 

Therefore at small the phase-space factor rapidly increases as 9i approaches 0° or 
180°. The functions Wt and Wtt become the dominant contributions in Eq. ([TI|) and 9i 
dependence of S is determined by xt and xtt- 

sin^ 9i + 2M^/P sm^ 9i 
Xt — — — , Xtt — ~ — • ("J^) 

2{smHi + M^/k^) 2{sm^ 9i + M^/k^) 

These are also rapidly changing functions, namely, xt increases from 1/2 to 1 and xtt 
decreases from 1/2 to when 9i approaches 0° or 180°. The product \l\^S determines the 
behavior of the the cross section shown in the top panels of Fig. (^). At large invariant 
masses the phase-space factor is independent of the angle. The 9i dependence of the cross 
section then comes only from the factors Xi in Eqs. ( p!2D and turns out to be rather fiat for 
this particular kinematics. 

Even if one moves away from the BH peak there is still a BH background. With increasing 
its relative importance increases as is seen from Figs. (^Q) . Therefore one has to choose 
a relatively large 9^ to obtain the situation in which the BH contribution can be ignored 
with respect to virtual CS. There is a relation between and the angle which we 
define as the photon scattering angle, where virtual CS is larger than BH by approximately 
one order of magnitude. The dependence is shown in Table | for two energies (for (almost) 
coplanar kinematics 9i is taken to be larger than 9bh to suppress the BH peak). As follows 
from Table | favorable conditions in the A region are at invariant masses up to ~ 200 MeV. 
At larger values of the BH background becomes comparable to or larger than virtual CS 
at all angles 9^. At an energy of 700 MeV larger invariant masses up to ~ 500 MeV become 
accessible. 
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B. Results for response functions in virtual Compton scattering 

In Figs. (^0) the response Ws = Wt + Wl and the ratios of RFs (or elements of the 
density matrix defined in Eqs. ([T8| , p!9|) ) are shown. The calculations are presented with the 
parameter set "B", as it gives better description of the photon asymmetry in real CS. 

First one notices from the angular distributions that at 9^ = 0° and 180° the interference 
terms pi_i and pio vanish, that can be explained from rotation symmetry around the beam 
axis. This property serves as an additional test of the model. 

In the A region (Fig. (^)) the transverse RF is dominant over the whole interval of M^. 
We studied the contribution of the additional coupling G^^k"^) in the 7A^A vertex in Eq. ( P5[ ) 
by varying within the limits ±g2- The effect turns out to be extremely small. This is 
due to small k'^ < 0.08 GeV^ which can be reached at this energy, and the dominance of the 
magnetic /^./m transition, which is very weakly dependent on G3{k'^). There is an effect of 
a exchange on the interference RFs, or on pi_i and pio (compare solid and dashed curves 
in Fig. (^). This is interesting in view of the very small contribution of sigma exchange 
to the cross section near the A peak. However, the region of above ~ 150 MeV will 
be difficult to access due to large BH contribution (see Table 1). The angular distribution 
of Ws at fixed is similar to real CS cross section da{'y)/dQ^, and is determined by the 
dominant Ml multipole and its interference with El and E2 multipoles. 

At a higher photon energy, shown in Figs. (^,^, the transverse response is still rather 
substantial. In Fig. (^) the effects of the important contributions are shown by switching 
off the corresponding process. The effect of the Di3(1520) resonance is clearly seen in a pro- 
nounced angular distribution while other contributions give rise to a flat background (dotted 
lines) with Qpio as an exception. This background in pi_i and ^pio is also independent of 
as is seen in the left panel. Other important contributions in this energy region are the 
5*11(1535) excitation and the a exchange. The 6*11 resonance contributes 10-20% to the cross 
section at = 300 MeV, and also to poo and 3fJpio at backward angles; its role increases 
with increasing invariant mass, a exchange, as a t-channel contribution, shows up at back- 
ward angles with the exception of S^pio, where it is seen at forward angles. There is also a 
small contribution to pi_i from the 1^33(1 700), however, its effect is not shown separately 
because the e.m. couplings of this resonance could not be fixed quite accurately from pion 
photoproduction. We should emphasize that all contributions add coherently in the total 
amplitude resulting in a strong interference. In particular, the P33(1232) resonance, through 
an interference, affects observables in the 2nd resonance region. 

For comparison we also present calculations with the parameter set "A" (long-dashed 
lines). In Ws large differences between the set "A" and set "B" appear at 6^ < 90° and 
large invariant masses, on average the differences are about 20%. The most striking effects 
can be observed in poo, Pi-i and 3?pio (compare solid and long-dashed hues in Fig. (|D). 
The big difference in pi_i at 90° is similar to that for the photon asymmetry in Fig. (^, 
while at larger angles the difference is reduced. As is seen from Fig. (|) for poo and 3?pio, 
the longitudinal response also proves to be rather sensitive to parameters 01^2 for the -D13. 

In principle one might expect a contribution of the Roper resonance, Pii(1440), to the 
longitudinal response. Data on pion electroproduction (see review [^) indicate however 
that the corresponding helicity amplitude S1/2 is very small, consistent with zero. The 
longitudinal coupling in the 'yNPn vertex has thus been neglected. The transverse coupling 
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is known to be small from real-photon data. 

In Fig. d^) the effect of G^{k'^) in the 7A^i5i3(1520) vertex is demonstrated. The dashed 
curves are calculated with = —3.3 which is chosen to reproduce the k"^ dependence of the 

1 /2 

ratio of the helicity amplitudes A1/2 and A3/2 (or the corresponding multipoles and 
M^L^) for pion electroproduction [^. As is seen, the G^^k"^) coupling has a considerable 



effect on Ws and thus on the cross section. The reason is that G^{k'^) strongly influences 
the electric dipole transition El and the latter in turn becomes the dominant term at large 
positive fc^. For the following consideration we choose y/s = {Mr is a resonance mass), 
then the 3-momentum squared of the virtual photon is 

k^ = kl -e = T^[(M^ + -e][{Mn- mf -k^]. (36) 

It is seen that k'^ decreases in the time-like region and vanishes at the kinematical limit, 
where fco = M^"^ = Mr — m. 

The behavior of RFs for e~e"*" production can be understood based on general properties 
of the resonance multipoles at small \k\ -C -R~^, where i? is a typical interaction radius. 
According to |^6| (ch. 6.2 and 6.3) one has for the e.m. transition with multipolarity j 



T{k- Mj) oc \k\^ , T{k; Ej) oc \kr'ko , T{k; Lj) -T{k; Ej) {j^f^' ■ (37) 

Therefore for the D13 resonance magnetic quadrupole (M2) vanishes, and electric dipole (£^1) 
and longitudinal dipole (LI) remain finite at fc^ = (M^""^)^. At the same time all multipoles 
for the P33(1232) resonance, magnetic dipole (Ml), electric and longitudinal quadrupoles 
{E2 and L2), are proportional to \k\ and go to zero. To have a more quantitative estimate 
we can use predictions of the non-relativistic quark model (e.g., Ref. |0). In Table |I| 



expressions are collected for the e.m. helicity amplitudes and the resonance pion-production 
multipoles. It follows from Table ^ that: 

i) Ai/2 and A3/2 for the D13 are finite at small \k\, while those for the P33 vanish; 

ii) if \k\ the ratio Ai/2/^3/2 for the D13 takes the value l/\/3 which coincides with 
the corresponding ratio for the P33 (in this and other SU{6) symmetrical models); 

iii) as /c^ increases towards the maximal value the spin-flip contribution (terms oc k'^) 
diminishes for the D13, correspondingly E2-. increases and M2- decreases (absolute values). 

Since the El transition dominates the D13 7*iV process, the coupling G^^k^) affects 
strongly the cross section in Fig. (|^). The transverse response Wt at the resonance position 
is roughly proportional to |-E1|^, while the longitudinal response Wl to |Llp|ii^lp. At 
small \k\ one can make use of Siegert's theorem (last relation in Eqs. (^7|)) which relates LI 
and El amplitudes. From this it can be seen that the ratio poo = Wl/{Wt + Wl) is equal 
to 1/3, irrespective of G^i^k"^). Other elements of the density matrix on Fig. (|^), except Q'pio, 
show a dependence on Gs{k^). 

It is of interest to compare Figs. (||J^) for the D13 with Fig. (^ for the P33 resonance. Due 
to a decrease of the Ml intensity as a function of M^ the response Ws for the A resonance 
falls off towards the maximal M^. The longitudinal response Wl contains the small resonance 
multipole L2 which also vanishes at \k\ 0, while the background receives nonvanishing 
contributions from LI multipoles (such as Lo+ and L2- in pion electroproduction). These 
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cause an increase of poo at large photon masses in Fig. (|) . This effect comes as a result of a 
balance between Gi{k'^) and G2{k'^) contributions in the 'fNA vertex in Eq. (|33D . Note also 
that in the quark model in Table |I| the electric quadrupole is absent, E2 = L2 = 0, which 
corresponds to a particular choice of the couphngs 2mgi = M^g2 = M^g^. 

The features of the e.m decay of the resonances in the time-like region are, in general, 
different from those studied in electron scattering, where |fc| is always larger than fco- For 
instance, in pion photo- and electroproduction interesting properties of the Z)i3(1520) reso- 
nance have been observed at —3 GeV^ < < (see, e.g. [0). In terms of the multipoles. 



the ratio E1/M2 —1 at large negative fc^, crosses zero at fc^ ~ —1 GeV^ and is about 
+2.1 at the real-photon point. It would be of interest to see if this ratio keeps on increasing 
at positive k"^ and becomes very large near A;^ = (M/j^g — m)^ = 0.3364 GeV^. 

Finally we note that the increase of Ws in Figs. above ^ 500 MeV is due to 
the fact that the e.m. form factors contain the p-meson propagator. As approaches 
the p-meson pole (which would be possible at photon energies above 1 GeV), the process 
7P — i> e~e~^p proceeds through creation of the p meson, as assumed in VMD models. 



IV. CONCLUSIONS 

We investigated virtual Compton scattering on the proton in the time-like region (7^ 
e~e~^p). When in an experiment the e~ and are not distinguished, the Bethe-Heitler- 
nuclear interference vanishes and the cross section is the incoherent sum of cross sections for 
the BH and the nuclear processes. We have shown that, contrary to common preconception, 
a considerable part of the phase space, mainly at backward angles for the virtual photon, 
is hardly contaminated by the Bethe-Heitler process. Under these conditions it is possible 
to decompose the exclusive cross section in terms of response functions. These are directly 
related to the photon density matrix pxy which characterizes polarization properties of the 
virtual photon. This offers a possibility to analyze e~e~^ production experiments also in 
terms of matrix elements pxy or, equivalently, tensor and vector polarizations of the time- 
like photon. 

Differential cross sections and matrix elements pw are calculated in a unitary K-matrix 
model which includes nucleon, mesons, and baryon resonances with masses up to 1.7 GeV. 
The model is tested for real-photon Compton scattering at energies up to 1 GeV. The agree- 



ment with data is reasonable and of a comparable quality to recent K-matrix calculation ||35 
and dispersion- relation analysis PT| . 

The dilepton production in the A-resonance region is dominated by the transverse re- 
sponse. At the largest at this energy certain elements pw show the effect of a exchange. 
The latter could be an approximation to various t-channel scalar- isoscalar exchanges (e.g., 
two-pion, e(760) and /o (400-900)). 

Photon energies corresponding to the 2nd resonance region, of about 700 MeV, allow one 
to explore higher photon masses up to ~ 500 MeV. Main mechanisms contributing here are 
the Di3(1520)-resonance excitation and, to a lesser extent, the 5*11(1535) resonance and a 
exchange. The transverse-transverse element pi„i strongly depends on the so-called off-shell 
parameters ai and 02 in the e.m. vertex of the D13. The longitudinal poo and longitudinal- 
transverse 3fJpio elements turn out to be sensitive to all contributions and with increasing 
photon invariant mass this sensitivity becomes more pronounced. At large the cross 
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section and the density matrix become strongly dependent on the part of the e.m. vertex of 
the 1^13(1520) which contributes only for virtual photons. 

Response functions are thus shown to be an important tool to distinguish between dif- 
ferent mechanisms in the e~e~^ production in the resonance region. If the resonance contri- 
butions can be separated, this will allow for studying their e.m. properties in the time-like 
region which may give an information supplementary to that obtained from the electron 
scattering. 
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APPENDIX A: RESPONSE FUNCTIONS IN HELICITY FORMALISM 

It is convenient to introduce the following set of polarization vectors for a time-like 
photon 

e*(0) = Jf^i\k\,0, 0, h) , e*(±l) = -^(0, ^1, ^, 0) (Al) 
which satisfy e*(A) ■ A; = and 

e*{X)-e{X') = -6,y, ^^'i^^VW = -9"" + ^ ■ (A2) 

A=0,±1 ^^-'7 

The RFs defined in Eq. ( 1T3|) can now be written as 

Wt = IY: \J- + + \J ■ e*{-l)\' = \ E |/+1A',AAP + |^1A',AAP , 

polar. ^ A,A'=±l/2, A=±l 

wl = \t. \j-^*m' = \ E i/oA'.AAr, 

polar. ^ A,A'=±l/2, A=±l 

Wtt = \ E 2^J-e%+l)[J-e*{-l)r = \ E 23? ,a/!ia',aa , 

polar. ^ A,A'=±l/2, A=±l 

Wlt=]T. 23? J ■ e*(0)[J ■ e*(+l) - J ■ e*(-l)]* 



4 , 

polar. 



- ^ 23?/o\, ;^^[/+iA/_AA - /-ia',aa] , (A3) 

A,A'=±l/2, A=±l 



where we introduced the helicity amplitude 
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(A4) 



which is a function of s,9^ and M^. The additional RF Wlt-, which defines the vector 
polarization Py = i-^tu of the virtual photon (see Eq. (p4D), can be calculated using the 
last formula in Eq. (|X3|), where the operation 3? ... should be replaced by .... 

Due to space-reflection invariance there are 8 (4) independent transverse (longitudinal) 
amplitudes in the sums in Eqs. ( [A 31 ). They can be chosen according to Table The remain- 
ing amplitudes are determined from the relation /_a'-a',-a-a = (—1)^ ~^ ~*-^~^^/a'A',aa [H 



In terms of fi of Table |T| one can rewrite the response functions as 



12 



Wr 



^ i=9 



2 

1 1 ; 



Wtt = '^{ fif; + f2f: + hf; + un ) , 

Wlt = 3? [/g (/i - /s) - /ro(/2 - U) + fuih - fr) - fUU - fs) 



(A5) 



For real photons due to the time reversal one has in addition = and /g = —/J, 
where the superscript "r" indicates the real CS limit: fyj^, xa^^^ ^7) = /a'A',aa('5, 0-y, = 0). 
In this case 



= -[ l/iT + 1/21^ + 21/3^2 + 21/ip + i/^f + 1/, 



r|2 ■ 



w^^T = 3^ [ /r(/[ + fD + /r(/2 - /e 



(A6) 



since /g, • ' ' > /12 vanish. Eq. agrees with definitions given in [2|,|C] up to a common 
normalization factor. 

In addition to the photon asymmetry in Eq. ( pQ]) we will also calculate the polarization 
of the recoil proton in real CS. It can be expressed through the helicity amplitudes as 



ppr 

y 



E 



4Wf A=±i/2,A,A'=±1 



9C> fr rr* 

iA'-l/2,AAiA'+l/2,AA 

1 



^>[/r(/[+/5)+/r(/6-/2) 



(A7) 
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TABLES 



^7 = 

My (MeV) 


320 MeV [M^ = 278 MeV) 
0" (deg) 


= 700 MeV (Al" 
(MeV) 6 


= 543 MeV) 

mm (deg) 


5 


30 


5 


30 


50 


90 


50 


60 


100 


120 


100 


70 


150 


150 


150 


90 


200 


170 


250 


110 






350 


130 






450 


150 






530 


170 



TABLE I. Minimal photon angle (defined in sect. Ill A) as function of photon invariant mass 



Di3(1520), 1=1/2 (proton) 



P33(1232), 1=3/2 



^1/2/^3/2 

EL (El) 
E{+ {E2) 
(Ml) 
Mi_ (Af2) 



2i{k^ - b^)aF{k^)ng 
-2iV3b'^aF{P)ng 

^(P-462)aF(P)/x,/ 



-2^\k\F{P)f,g 
1 

V3 



TABLE IL Electromagnetic helicity amplitudes and pion-production multipoles for Di3(1520) 
and ^33(1232) resonances in model [47]. F{}?) = (1 + a^P)~^, 5^ = rnqk^/gq, and mg,fiq,gq and 
a are respectively constituent quark mass, scale magnetic moment, gyromagnetic ratio, and scale 
parameter. The constant / = {it /ko)^/"^ {k^ymTT^N / ATrgTrMjiT^^^y^'^ relates helicity amplitudes to 
e.m. multipoles, where k^ is equivalent real-photon energy, is pion 3-momentum, F^at {Ttot) is 
one-pion (total) decay width of the resonance. 



/l 


f+1 


fl/2,- 


n- 


fl/2 


/5 


/-l-l -1/2, +1-1/2 


/9 


/o- 


f 1/2, +1+1/2 


/2 


f-1 


-1/2,- 


n- 


fl/2 


/6 


/-I +1/2, +1-1/2 


/lO 


/o- 


-1/2, +1+1/2 


/3 


f-1 


-fl/2,- 


n- 


fl/2 


/7 


/-I -1/2, +1-1/2 


/ll 


/o- 


-1/2, +1-1/2 


h 


Ui 


-1/2,- 


n- 


fl/2 


/8 


/+1 +1/2, +1-1/2 


/l2 


/o- 


fl/2, +1-1/2 



TABLE in. Lidependent helicity amplitudes in the 7p — > e e~^p reaction. 
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FIGURES 




FIG. 1. Tree-level graphs describing the 7A'^ — > e~e^N process. Diagrams a f correspond 
to virtual CS on the nucleon, g and h to the Bethe-Heitler contribution. Single (double) solid 
line depicts a nucleon (baryon resonance), wavy line a photon, dashed line tt*^, a and rj mesons. 
Diagram f describes a possible contact ^^*NN vertex. 
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FIG. 2. Differential cross section (top), photon asymmetry (middle), and recoil-proton polar- 
ization (bottom) for real-photon Compton scattering as function of photon energy. Solid (dashed) 
lines are the present model calculations with the parameter set "A" ("B"), dotted lines are calcu- 
lations from 1 31]. Data for cross sections are: left panel - • Mainz ||3^, ★ Cornell |3^, A MIT ]37| , 
o Bonn [||, o Tokyo [||; right panel - • Bonn ||, o Tokyo □ Moscow (Lebedev Inst.) gl|. 
Data for photon asymmetry: o LEGS Collaboration |42|, • Erevan [43|; and for proton polarization: 
o Tokyo Q. 
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FIG. 3. Exclusive differential cross section in the cm. as function of lepton "sharing" angle 9i 
in coplanar kinematics. Left panel corresponds to E!:^^ = 320 MeV, right panel to EI^^ = 700 MeV. 
The virtual-photon angle is 6j = 135°. Dashed and solid lines are respectively the Bethe-Heitler 
and the virtual Compton scattering (parameter set "A") contributions. 
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FIG. 4. Exclusive cross section at an incoming energy of 700 MeV in non-coplanar kinematics. 
Notation for the cm'ves is the same as in Fig. @. 
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FIG. 5. Response Ws = Wt + Wl (top) and polarization density matrix versus virtual-photon 
invariant mass My (left) and photon angle 6^ (right). Incoming photon energy is 320 MeV. Solid 
lines are calculated with the total amplitude, dashed lines without a exchange. 
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FIG. 6. The same as in Fig. @ but at photon energy 700 MeV. Sohd hnes show the calcu- 
lation with the total amplitude (parameter set "B"), dotted lines without £)i3(1520) resonance, 
dash-dotted lines without 5ii(1535) contribution, short-dashed lines without a exchange, and 
long-dashed lines are calculated with the total amplitude (parameter set "A"). 
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FIG. 7. Response Ws = Wt + Wl (top) and pw calculated for different values of in the 
^NDi^ vertex Eq. (p3|). Solid lines: = 0, dashed lines: 33 = —3.3 (from pion electroproduction), 
and dotted lines: 53 = —7.9. Couplings gi and g2 are respectively 7.25 and 7.9. 
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